Statistical physics of the ScheUing model of segregation 
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We investigate the static and dynamic properties of a celebrated model of social segregation, 
providing a complete explanation of the mechanisms leading to segregation both in one- and two- 
dimensional systems. Standard statistical physics methods shed light on the rich phenomenology of 
this simple model, exhibiting static phase transitions typical of kinetic constrained models, nontrivial 
coarsening like in driven-particle systems and percolation-related phenomena. 



Individuals with similar ideas, habits or preferences 
tend to create cliques and cluster in communities, which 
results in segregation at the scale of the society. A first 
theoretical analysis of segregation phenomena in social 
environments was performed by T.C. ScheUing, who de- 
fined a model in which agents divided in two species move 
on a checkerboard according to a given utility function 
Within this simple model, ScheUing was able to show 
that segregation occurs even when individuals have a very 
mild preference for neighbors of their own type, as long 
as they are allowed to move in order to satisfy their pref- 
erence. 

Interest for Schelling's results has grown recently among 
social scientists 0, Ijl, mathematicians J3| and statisti- 
cal physicists [1]. In particular, Ref. [J suggested re- 
lations with models of binary mixtures in physics. In- 
deed, the way in which segregation takes place in the 
ScheUing model resembles the coarsening processes gov- 
erning phase separation kinetics ^, 0]- However, apart 
from a qualitative picture, Ref. [2] does not provide any 
quantitative measure of coarsening from which the scal- 
ing behavior of the segregation process could be inferred. 
On the other hand, looking at the non-equilibrium pro- 
cess of such a model, several questions spontaneously 
arise: Does the system fully segregate? What are the 
properties of the stationary state? If coarsening takes 
place in this system, does it fall in one of the known 
universality classes? 

In this Letter, we shed light on the static and dynam- 
ical properties of the ScheUing model in one and two di- 
mensions, applying methods taken from different fields of 
statistical physics, such as (vacancy mediated) coarsen- 
ing dynamics J], diffusion-annihilation particle systems 
0, and kinetic constrained systems [l^. 
The model is the same investigated in Ref. individu- 
als are initially distributed at random on a line {d = 1) or 
a square lattice (d = 2) of iV = L'' sites. The occupation 
of each site i is described by a spin variable cr^, taking 
values (Ti = if the site is empty and = ±1 if the site 
is occupied by an individual of type ±1. Let po — Nq/N 
and p± = N±/N be the densities of vacancies and of oc- 
cupied sites, of either type, respectively. 
The general principle governing the dynamics is that an 
individual can tolerate at most a fraction f — f* of his 



(Moore) neighbors to be different from him. li f < f * 
the individual is happy or has utility 1 otherwise he feels 
unhappy (utility 0). Unless otherwise stated we consider 
/* = 1/2, as in the original ScheUing model. Two types 
of dynamics are possible: constrained dynamics (called 
"solid model" in Ref. Q), where only unhappy individ- 
uals are allowed to move, as long as they are able to find 
a vacancy where they can be happy; and unconstrained 
dynamics (called "liquid model" in ;2i|) where agents are 
allowed to move to vacancies as long as their situation 
does not get worse off. In both cases, following Ref. 0], 
we assume an infinite range dynamics, i.e. individuals 
can move to any suitable vacancy, irrespective of dis- 
tance. This spin-vacancy exchange dynamics conserves 
only globally the magnetization, but not locally. Natu- 
ral quantities to characterize the state of the system are 
the densities of unhappy sites u{t) and of interfaces n(i), 
defined as the fraction of neighboring spins of opposite 
type. 

We stress that, while it might be appealing to introduce 
a notion of energy, as e.g. the number of unhappy indi- 
viduals, this might be confusing as the dynamics is not 
derived from an energy functional, as in physics. Individ- 
uals move solely in order to maximize their utility, with 
no regard to the welfare of the fellow neighbors. For ex- 
ample, even if a displacement is beneficial to the mover, 
this might make some of his neighbors unhappy with the 
composition of their new neighborhood. So a move may 
cause an increase of the number of unhappy individuals. 
In the following, we focus on the static properties of the 
constrained model and the dynamical properties of the 
unconstrained one, showing that statistical physics allows 
one to understand many aspects of the rich phenomenol- 
ogy of the model. We shall use the Id case to uncover 
the main properties of the segregation process and show 
that this provides key insight on the behavior of the 2d 
case. 

The constrained model - Let us first discuss the case 
where only "condensation" moves, which increase the 
utility of the moving agent, are allowed. At each time 
step, an individual is drawn at random and swapped 
with a randomly chosen vacancy if his utility increases. 
If there is no such vacancy, the agent is not displaced. 
This dynamics converges in a finite time to frozen states. 
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FIG. 1: a) Density of unhappy spins Uoa in the blocked con- 
figurations of the constrained Id model as a function of the 
density of vacancies po and for different values of magnetiza- 
tion m > 0. Solid lines are obtained with the transfer matrix 
approach for the case of equal a priori weight; the dashed 
lines are the results of the dynamical computation, b) The 
quantity Uoo(po) computed dynamically for a 2d system of 
size L = 128. 



called (Myopic) Nash Equilibria (NE), where no agent 
can find a better location in terms of utility given the 
location chosen by the others [1, |3| • Segregation occurs 
up to a finite length, both in d = 1 and 2. 
In Id, it is possible to investigate the blocked configu- 
rations by means of a static approach [13] , that consists 
in considering the ensemble of such configurations with 
equal a priori weight. The partition function for this en- 
semble is 



ZL{h,i2,P) 



/iAf((j)+^7V(<T)+/3(7(o-) 



(1) 



where C is the set of all blocked configurations a = 
(cTi, . . . , (tl) of size L, with at — 0, ± (with periodic 
boundary conditions), M{a) = '^i'^i is the magnetiza- 
tion, N{a) = J2i ^1 is the number of individuals and 
U{a) is the number of unhappy individuals. The coef- 
ficient of the term proportional to v^'^p+~p-^w'"^''+^p-^ 
in Zi^ilnv, Inu", 0) is the number of NE for a chain with 
Lp+ individuals of one type and Lp_ individuals of the 
other type. Zl can be computed using the transfer ma- 
trix technique, following the same approach of Ref. 11 1. 
In brief, one can write Zl = TrT^, where the transfer 
matrix T relates the statistical weights of configurations 
of length £ -|- 1 to those of configurations of length I. In 
order to build configurations without constraints, T must 
have elements T^<^i-i-'^i)^i<^t''^w) which depend on three 
consecutive spins, i.e. T is a 9 x 9 matrix 14 1. 
It is important to realize that NE can be of three types 
only. Either a configuration is blocked because there are 
no unhappy sites, or because there are unhappy sites 
with, say, ai = + but there are no vacancies which could 



accommodate them (such as e.g. . . . -(- -I- . . .). This 
implies that the partition function has the form 

(/i, Ai, /3) = Tr + Tr + Tr Ti' (2) 

where Tq generates all configurations with no unhappy 
sites, while 7^ and those with unhappy sites of type 
+ or — , respectively, but no corresponding vacancies. 
For L ^ 1, each term in Eq. ([2]) is dominated by the 
largest eigenvalue \q{h, ^, fi) of the corresponding matrix 
Tq (q = 0,±). Hence logZL/L = max^ log Aq(/i, /i, /?). 
The number of NE with m = p-\- — p- and a density po 
of vacancies is given by e^^'^"^'P°''^^'^^ where the entropy 
S{m,po;l3) = max,,,^,, [log A,(/i, /u, /?) - hm ~ p{l - po)] 
is obtained from Ag via Legendre transform. This allows 
us to access the statistical properties of NE depending 
on the various parameters of the problem. 
In particular, it is interesting to look at the behavior of 
the static density Uoc — dpS{l3 = 0) of unhappy individ- 
uals as a function of po- For po > Po("^) solution 
is dominated by the configurations with no unhappy site 
(uoo = 0, Ao > A±) but the converse is true for small 
density of vacancies [po < Po('^)]- Indeed the density 
Woo(po) of unhappy sites features a jump at pp (Fig.[T^). 
Such first order phase transition is typical of cases where 
the partition sum can be separated in different compo- 
nents as in Eq. For m > the leading contribution 
to Zl comes from the term 7^, implying that NE are 
characterized by a finite density of majority type indi- 
viduals who cannot find a suitable location. For m — Q, 
the symmetry between the two types of agents is spon- 
taneously broken for po < Po(0), in the sense that a ran- 
domly drawn NE will typically have a fraction Uoo > 
of unhappy individuals of one type, all individuals of the 
other type being happy. 

The results of the static approach qualitatively agree 
with simulations, where blocked configurations are se- 
lected by the dynamics, starting from a random initial 
condition (see Fig. [Iji). When we start from m = 
p+ — p_ =/: 0, we again observe a non-zero fraction of 
unhappy individuals (uqo > 0) of the majority type, for 
small enough po. At odds with the static computation, 
i) the transition is continuous rather than discontinuous 
and a) when m = no spontaneous symmetry breaking 
is observed. [3] [Hi 

The same qualitative picture applies to the original two 
dimensional case, with Moore neighborhood. In d = 2, 
the transfer matrix approach to the static problem be- 
comes unfeasible. Progress can be done by MonteCarlo 
simulations of an appropriately defined auxiliary model 
[l^. This shows the presence of a discontinuous transi- 
tion in Uoo as a function of po, similar to the one occur- 
ring in Id (not shown). Moreover, Fig.[T|3 shows that the 
Moo(po) curves produced by the dynamics of Schelling's 
model have the same qualitative behavior of the Id case. 

In the unconstrained model, at each time step a ran- 
domly chosen agent is relocated into a randomly chosen 
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FIG. 2: Scaling behavior of the density of interfaces n{t) ob- 
served in the hmit of low density of vacancies po for a Id 
system of size L = 10® and zero magnetization m — 0. The 
inset shows the behavior of the stationary value of the densi- 
ties of unhappy spins and interfaces as a function of po- 



vacancy (long-range diffusion) when his utility does not 
decrease. No blocked configurations exist since individ- 
uals can always be relocated. For this reason, the sys- 
tem can enter a stationary regime in which the densities 
of unhappy individuals u{t) remain constant on average. 
We are thus interested in understanding the asymptotic 
properties of the dynamics, and characterizing the evo- 
lution to the stationary state. We study these properties 
in ci = 1 , 2 as a function of the density of vacancies po , 
assuming m = 0. 

In d = 1, for any pq, u{t) has an initial decay followed 
by a stationary state with asymptotic value Moo(po)- The 
initial decay is exponential for large densities of vacancies 
{Pq — > 1), while it becomes power-law in the limit of van- 
ishing Pq. The density n(t) of interfaces between domains 
of unlike spins exhibits a similar behavior. Fig. [2] shows 
that for Po ^ 1 the interface density follows the scaling 
law n{t) ~ poi^ipot^^^) with ^{x) ^ for a; ^ and 
ipix) const for a: — > cxd. A similar power-law behavior 
{u{t) ~ t~^^'^) is found for u{t) in this limit, but u{po,t) 
does not satisfy any scaling law. 

The behavior of the stationary densities as functions of 
Po is reported in the inset of Fig. [21 Both Uoo and rtoo 
have a maximum at some value < pg < 1, then van- 
ish in the limits po 0,1. In particular, for po 1, 
Moo(po) vanishes as (1 — po)^ , while in the opposite limit 
n-oo Po and Moo ~ Po- 

An explanation for the observed behavior of for 
Po ~ 1 can be given in terms of a phenomenological dy- 
namics for u{t). This takes the form u = —au + b where 
a and b describe processes which annihilate or create un- 
happy sites, respectively. For po « 1 an unhappy indi- 
vidual with high probability becomes happy, once it is 
displaced, hence a ~ 1. Unhappy individuals are created 

in processes such as. ..-I s-...-) — . . . where a 

displaced site leaves one of his neighbors in an uncomfort- 



able neighborhood. For po — l, the leading contribution 
is given by initial configurations of three occupied sites, 
whose probabifity is 4[(1 — po)/2]'^, where the factor 4 
accounts for the degeneracy of the possible types. Hence 
u(t) decays exponentially to Uoo = b/a ^ (1 — po)^/2, 
which agrees very well with numerical simulations. 

A more elaborate analysis is instead necessary in or- 
der to explain the observed phenomenology in the limit 
Po — > 0. It is convenient to think in terms of interfaces 
between clusters of opposite spin values. These can be 
thought of as particles, like in diffusion- limited annihi- 
lation processes. For two adjacent clusters of opposite 
spins, the leading process is the diffusion of interfaces, 

e.g. -I-H >-h + >-h + H . The other 

relevant processes are the creation and annihilation of 
clusters of size one, i.e. unhappy spins. The annihilation 
rate of unhappy spins is obviously cx n^, whereas the cre- 
ation of interfaces involves "vacancy mediated" processes 

of the type -I-H > + — 0— . Here the displacement 

of an individual makes one of his neighbors unhappy, and 
the latter can in turn move in the bulk of a cluster of indi- 
viduals of the other species. In the stationary state, the 
creation rate is proportional to the probability to find 
an empty site close to an interface, which is pon. The 
balance condition pon « implies that the stationary 
interface density is Uoo oc po and hence that the creation 
rate is proportional to pg. 

The solution to a problem of diffusion-limited annihi- 
lation with input of particle pairs is well known [9] and it 
leads to the following scaling form for the density of in- 
terfaces, n{po,t) ~ pQ^^tl){pot^^^), with tpix) — > const 
for a; — > oo and il^{x) — + x^^^^ for x ^ 0. The in- 
jection rate pg implies 5 — 1 and A = 2, in perfect 
agreement with numerical simulations. At po = 0"*" (sin- 
gle vacancy limit), the dynamics is purely diffusive, thus 
?i(po = 0+,t) ^ with z ~ 5 A ~ 2. In other words, 

the Po — > limit can be considered as a critical point for 
the one-dimensional unconstrained Schelling model. In 
this respect, the Schelling model for pg > produces a 
dynamics with the same universal features of finite tem- 
perature coarsening dynamics in Ising like models with 
non-conserved order parameter . 

Contrary to n{po,t), the density of unhappy individu- 
als does not follow a scaling form. Indeed one can write 
u{po,t) = —dn{po — 0,t)/dt + cpgn, with c > a con- 
stant. The first term is the probability to have a cluster 
of size one in the bulk of larger domains, whereas the sec- 
ond arises from the creation of unhappy individuals close 
to interfaces (see above). For early times (t <^ pg ^, as 
suggested by dimensional analysis) the first term domi- 
nates, hence u{po,t) ^ t^^^^, whereas u(pg,t) converges 
to nl^ipo) for t > po^ 

We now turn to discuss the two-dimensional uncon- 
strained dynamics, that was already studied in Ref. [1], 
where it was suggested that some coarsening phe- 
nomenon should take place. We remark that, in the dy- 
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FIG. 3: Density of interfaces n{t) in a 2d system of size 
L — 10"^: a) /* — 1/2, the coarsening process is present for 
Po < Po; b) the systems does not order for some values of the 
threshold /* / 1/2. 



state would be optimal, but it cannot be reached dy- 
namically, as the system remains trapped in a disordered 
state. 

Conclusions - In summary, we have investigated both 
analytically and numerically the static and dynamic 
properties of the Schelling model of segregation in one- 
and two-dimensional systems. The constrained version 
of the model presents non-trivial static properties char- 
acterized by the existence of a transition with symmetry 
breaking, whereas the unconstrained dynamics exhibit 
coarsening typical of systems with non-conserved order 
parameter. A further sharp transition takes place as the 
tolerance threshold (/*) of individuals gets either very 
large or very small, with the system being trapped in 
disordered dynamical states. 



namics, voids are expected to be randomly distributed 
in the system. If po is too large, a cluster of void sites 
starts percolating throughout the system, preventing the 
growth of the clusters of spins ^ beyond a finite size. In 
a d = 2 square lattice with Moore neighborhood, the per- 
colation transition takes place at po ~ 0.407 [5|. Since 
here voids can move we expect this value to be only a 
lower bound for the real transition, that takes place at 
some 0.45 < po < 0.5. For po > po we find that clusters 
grow up to a finite size de pen ding on po but indepen- 
dent of the system size A'^ [13[. In contrast, below the 
void percolation transition, we observe convergence to 
a (quasi) ordered state, with two domains spanning the 
whole system. A clear coarsening process with the typ- 
ical length of clusters growing as t^^^ with z = 2 (see 
Figure [3^) is observed for small po. Even if the magne- 
tization is globally conserved, the value z = 2 correctly 
describes a coarsening process with non-conserved order 
parameter in agreement with our findings in Id and with 
renormalization group results for domain-growth scaling 
in the presence of long-range diffusion Q . 

Finally, it is interesting to discuss the behavior of the 
unconstrained model as a function of the maximal frac- 
tion /* of unlike neighbors which individuals tolerate (so 
far equal to 1/2) in rf > 1. Figure [Sb shows that, for 
Po ^ 0, the coarsening process is robust against the vari- 
ation of /* in a wide range around 1/2, as also found in 
Ref. 0. Segregation takes place even if individuals are 
satisfied with as few as 2 neighbors of their own type. If 
instead individuals are as tolerant as to be satisfied with 
just one like neighbors, then no coarsening and no segre- 
gation takes place and the system remains in a dynamic 
disordered state. No coarsening takes place also when in- 
dividuals are extremely intolerant (see Fig. [SId). This is 
somewhat remarkable, because indeed a fully segregated 
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